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Abstract 



0^ ' In this article we continue our investigations of one particle quantum scattering 

^^ ■ theory for Schrodinger operators on a set of connected (idealized one-dimensional) 

wires forming a graph with an arbitrary number of open ends. The Hamiltonian is 
Mh. given as minus the Laplace operator with suitable linear boundary conditions at the 

vertices (the local Kirchhoff law) . In "Kirchhoff's rule for quantum wires" [J. Phys. 
A: Math. Gen. 32, 595 - 630 (1999)] we provided an explicit algebraic expression for 
r^ ' the resulting (on-shell) S-matrix in terms of the boundary conditions and the lengths 

^^. of the internal lines and we also proved its unitarity. Here we address the inverse 

problem in the simplest context with one vertex only but with an arbitrary number 

of open ends. We provide an explicit formula for the boundary conditions in terms 

^^ ' of the S-matrix at a fixed, prescribed energy. We show that any unitary nxn matrix 

may be realized as the S-matrix at a given energy by choosing appropriate (unique) 
boundary conditions. This might possibly be used for the design of elementary gates 
in quantum computing. As an illustration we calculate the boundary conditions 
associated to the unitary operators of some elementary gates for quantum computers 
and raise the issue whether in general the unitary operators associated to quantum 
gates should rather be viewed as scattering operators instead of time evolution 
operators for a given time associated to a quantum mechanical Hamiltonian. We 
also suggest an approach by which the S-matrix in our context may be obtained from 
"scattering experiments" , another aspect of the inverse problem. Finally we extend 
our previous discussion, how our approach is related to von Neumann's theory of 
selfadjoint extensions. 
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1 Introduction 

In this article we continue our investigation of quantum mechanical scattering theory on 
a set of connected wires idealized to have zero diameter and with an arbitrary number of 
open ends [gO| . The dynamics is given by the one-dimensional Laplacian with arbitrary 
boundary conditions at the vertices, which makes it a selfadjoint operator. We view this 
as an idealized version of a thin conducting wire with electrons in the conductance band, 
i.e. we have set h = 2m* = 1 with m* being the effective mass and we have neglected the 
spin. Also if the wires are sufficiently thin, transversal modes may be neglected and only 
longitudinal modes, i.e. motion in the direction of the graph, should play a role. Below 
we will comment briefly on the associated matheinatical problem. A concrete example 
we have in mind are nanotubes, which are considered to be interesting candidates for the 



design of microscopic electronic devices (see e.g. |52, 35, SC] for recent, possibly relevant 
experiments). In fact, the possibility of connecting nanotubes of different diameter and 
chirality has created considerable interest recently (see e.g. |17, ^ |^, ^, ^, ^). Alter- 



natively such graphs could be realized in terms of grooves etched on suitable surfaces 
or in terms of strings of atoms forming a (planar) graph deposited again on a suitable 
surface [[|, g, |0|, see also |l], |3|]. 



For our graph model we provided an explicit algebraic expression for the S-matrix at 

a given energy in terms of the boundary conditions and the lengths of the internal wires. 

We proved unitarity so this provides an example for the Landauer-Biittiker formalism 

M|, p, see also e.g. |10] for an extensive discussion. We note that the transmission rate 



for the connection of two nanotubes with different diameter has been calculated in |51]. 
The influence of disorder on the conductance characteristics of nanotubes was studied 
ini. 

In this article we will also address the inverse problem, namely the determination of 
the boundary conditions given the S-matrix at a fixed energy. We will only cover the case 
of one vertex only, since the general case with internal wires and several vertices is much 
more difficult. Note that another kind of inverse problems was considered previously in 



0. 



We will use this discussion of the one vertex case to indicate possible applications 



in the design of quantum computers (see e.g. |g^, ^, 48 1 and the literature quoted 
there). Recall that in classical network theory one uses the theory of unitary matrices 
to describe (deterministic) input-output relations and one even speaks of scattering 
matrices (see e.g. [^, ^]). For quantum computers we advocate the idea that one 
should use the theory of quantum scattering and the associated notion of a unitary S- 
matrix to formulate (probabilistic) quantum input-output relations. This differs from 
standard discussions on this subject in which the unitary matrix is considered to be 
the unitary time evolution operator for a given, fixed time and with a Hamiltonian as 
infinitesimal generator, which describes the dynamics. Also we consider the rule by 
which the connection of gates is supposed to correspond to matrix multiplication of the 
associated unitary matrices not to be quite convincing. In fact, connecting two gates 
means that one has a coupled system for which the relevance of the Hamiltonians of the 
two subsystems remains unclear. So from the point of view of information transmission, 
we consider it at least natural to ask for the relevance of scattering matrices in the 
context of quantum computation. In fact, it is our understanding that most of the 
present experimental designs for quantum gates describe scattering experiments. In the 
present context of quantum wires the incoming signals would be plane waves at a fixed 
energy in each of the wires and which will be scattered at the vertex into outgoing plane 
waves again in each of the wires and of the same energy. 



In this context it is worthwhile to mention the connection between the time depen- 
dent quantum mechanical description and scattering theory. The latter describes the 
long time behaviour of quantum evolution. Thus the use of scattering matrices instead 
of unitary time evolution operators can only be adequate provided the "tact frequency" 
of the quantum computer is not too high and this depends on its physical realization. 
However, a discussion of this issue as well as the question of how to incorporate this 
proposed realization of quantum gates into the general context of quantum computation 
with its notion of entanglement lies outside the scope of the present article. 

The discussion below will be based on our result that in the situation of a single 
vertex quantum wire any unitary matrix may be uniquely realized as the S-matrix at a 
given energy by choosing appropriate boundary conditions. In spirit this is similar to 
the discussion in e.g. |^3 by which any unitary matrix in principle can be realized with 
suitable beam splitters and phase shifters. As an illustration we will calculate these 
boundary conditions for some simple gates. We leave out the question of how to imple- 
ment these boundary conditions at the vertices of concrete physical wires. However, in 
the present context we will discuss the general experimental difficulty associated with 
designing quantum gates with a prescribed unitary matrix. It has to do with the well 
known fact that usually only the absolute values of matrix elements of the S-matrix are 
observable (conventionally given in terms of cross sections). We will suggest an "exper- 
imental" prescription, how one may also determine the phases of the S-matrix elements 
by suitable experiments. This method may also be applied to potential scattering theory 
for Schrodinger operators p^]. 



We conclude this article by elaborating in more detail than in |20| on the relation of 



our discussion of selfadjoint extensions with the corresponding theory of von Neumann. 

2 The inverse problem for a quantum wire with a single 
vertex 

In this section we first discuss the inverse problem for n (infinite) wires ending at a 
single vertex, i.e. we will determine the boundary conditions given the S-matrix. We 
briefly recall the set-up for this situation. Let the Hilbert space be given as 

w = ©r=iWi = ©r=i^^([o,oo)). 

Intuitively the n origins are to be identified. This will be achieved in a moment 
by describing the quantum mechanical dynamics. Elements ip ^ 7i will be written as 
V' = (V'l, '02, ■■■,'4'n) and derivatives as ^p' = {tp[,ip2, •••, V'n) We will call ipj the component 
of ip in channel j. The scalar product in TC is 

n 
i=\ 

with the standard scalar product on L^([0,oo)) on the right hand side. Now the dy- 
namics will be given by the Laplace operator away from the origin supplemented by 
suitable boundary conditions at the origin. For this we proceed as follows. Consider the 
symmetric operator A*^ on "H, such that 



with domain of definition P(A*^) being the set of all V with ipjip'^tp" £ TC and such 
that ipi{0) = ip'iiO) = for all 1 < i < n. All selfadjoint extensions are now given as 



follows. Let T> C TC he the set of all ip such that ip , ip' , ip" are all in Ti.. On T> consider 
the following skew-Hermitian quadratic form given as 



with the Laplace operator A = —cP/dx'^ considered as a differential operator. Obviously 
O vanishes identically on T>{A^). Any self-adjoint extension of A*^ is now given in terms 
of a maximal isotropic subspace of P, i.e. a maximal (linear) subspace on which fi 
vanishes identically. To find these maximal isotropic subspaces, let [] : V ^ C^" be the 
surjective linear map which associates to ip and ip' their boundary values at the origin: 

[V] = (Vi(o),...^„(o),^;(o),...<(o)f = ( ^/°)^ 

Here T denotes the transpose, so [V'], ^(0) and ^'(0) are considered to be column vectors 
of length 2n and n respectively. The kernel of the map [] is obviously equal to I?(A'^). 
Then we have 

where ( , )c2n now denotes the scalar product on C^" and where the 2n x 2n matrix J 
is the canonical symplectic matrix on C^": 

Here and in what follows I is the unit matrix for the given context. 

To find all maximal isotropic subspaces in D with respect to fi it therefore suffices to 
find all maximal isotropic subspaces in C^" with respect to co and to take their preimage 
under the map [ ]. Let the linear subspace Ai = MiA, B) of C^" be given as the set of 
all [tp] in C^" satisfying 

AiPiO) + BiP'iO) = 0, (2) 

where A and B are two n x n matrices. If the n x 2n matrix (A, B) has maximal 
rank equal to n and if AB^ is selfadjoint then M.{A,B) is maximal isotropic and in 
this way one obtains all maximal isotropic subspaces in C^". Two maximal isotropic 
subspaces M{A'^B') and M{A,B) are equal iff there is an invertible n x n matrix 
C such that {A\B') = (CA,CB). In this sense we will say that the n x 2n matrix 
{A, B) is uniquely fixed by the boundary condition. Let T){A^ B) denote the inverse 
image of M{A,B) under the map []. Then P(A°) C 'D{A,B) C P and A f T>{A,B) 
is the core for a selfadjoint extension of A^ and which is denoted by A{A,B). This 
one-dimensional set-up could be a limiting case of a more realistic set-up with wires 
having non zero diameter. More precisely consider a graph to be a fixed subset of M^ 
such that in particular the angles between the edges ending at any vertex are fixed. 
Now form a tubular e-neighborhood of this graph, i.e. a real system of connected wires 
of diameter e. In mathematical terms this means that the graph is a deformation 
retract of any of these tubular neighborhoods. Consider the Laplace operator on this 
3-dimensional set with Neumann boundary conditions its boundary, i.e. on its surface. 
Now the mathematical question arises whether for e J, the resolvent (i.e. the Green's 
function) converges suitably to the resolvent of a Laplace operator on the graph and 
which is of the above type. If this is the case its boundary conditions must then be 
given in terms of geometrical data, i.e. on how this graph lies in M^. In other words 



and as already mentioned in the introduction, when e | the transversal modes (except 
0) should tend to infinity while the "longitudinal" modes should converge to those of 
a suitable Laplace operator on the graph. For some simple cases this is indeed known 
to be true (see for example |Q). The S-matrix S{E) = Sa,b{E) for any fixed energy 
ii^ > is now given as follows. It is an n x n matrix whose elements are defined by the 
following relations. Consider the plane wave solutions ip {-^E), \ < k < n oi the form 

V,'(x, E) = 6jke-'^^ + S,k(.E)e'^- (3) 

and which satisfy the boundary conditions. Then the S-matrix at energy E for the 
boundary condition {A, B) is given as 

Sa,b{E) = -{A + iVEB)-\A-iVEB) 

= -{A^ -^^/EB^){AA^ +EBB^)-^{A-iVEB) ^' 



and is unitary and real analytic in £' > 0. Note that Sa,b{E)~^ = Sa,-b{E). In ||2^ 
it was shown that the knowledge of the S-matrix at a fixed energy, Eq say, uniquely 
fixes the boundary conditions. Here we will provide an explicit solution to this inverse 
problem. For known S = Sa,b{Eo) let 

^' = -^(^-1), B' = -^{S + l). (5) 

z li\J Eq 

Then by (|) one has (A', B') = {CA, CB) with C = {A^ i^^B)-^ . This proves that 
S = Sa,b{Eq) indeed uniquely fixes the boundary condition. We can even say more. 
Indeed, observe that {A',B') defined by ^ satisfies the two relations 

A'B'^ = — \= (S^ -S), A' + iJWqB' = I, 

where the first relation is a simple consequence of unitarity. From the first relation we 
can also directly read off that A'B'^ is selfadjoint and from the second that [A', B') has 
maximal rank equal to n. Thus we have proved the following 

Theorem 1. For any unitary n x n matrix S and any energy Eq > there are unique 
boundary conditions {A,B) such that SA,BiEo) = S. More generally the relation 



Sa,b{E) = ({^ - ^/E'o)S + {Ve + ./Wo) , 
■ i^i^/E + /E^)S +iVE- /Eo)) 

holds for all E > 0. 

Indeed, the relation (|6|) follows from Sa,b{E) = Sa',b'{E) for all i? > and from 
relation (0). Since Eq was arbitrary this in particular gives the relation 



SaME) = {{Ve - /E^)SaMEo) + {^ + VE^)) 
[Ve + /E^)Sa,b{Eo) + {Ve- Ve^)^ 



for all E and Eq > 0. In particular this relation may be used as a test to verify if the 
S-matrix indeed results from boundary conditions in the way described here. 

As an illustration with possible applications to quantum computers in mind we 
consider the some elementary gates usually discussed in this context (see e.g. M). Thus 



we will view a vertex with n wires entering as a quantum gate with the wires viewed 
as channels. In particular the output channels are the same as the input channels. 
Therefore in addition to quantum transmission from one channel into a different one 
there will also be reflection from a channel into itself. These reflection amplitudes 
correspond to the notion of backpropagation in classical information theory. 

It is important to realize that here we are dealing with a one-particle theory. In fact, 
the notion of a particle is introduced here for the following reason. In the concrete setup 
of classical computers one deals with signals, which are localized in space and time. 
But in quantum (field) theory and with the appropriate modifications this is one of the 
notions used to describe particles in terms of wave packets. In the quantum mechanical 
formulation incoming signals at gates result in outgoing signals and this corresponds to 
looking at scattering theory and the associated unitary S-matrix. In particular in our 
approach the vertices act like classical external potentials for which the associated quan- 
tum scattering theory has been discussed extensively (see e.g. |4^). Here the dimension 
of the space of incoming plane waves at energy E is equal to n, the number of wires 
entering the vertex. This contrasts with the usual discussion of quantum gates with n 
input and n (different) output channels, where the Hilbert space has dimension equal to 
2" and replaces the set of classical information containing n bits. This dimension reflects 
the tensor product structure (n factors, each being two-dimensional) of the underlying 
theory. The tensor product structure is essential for the notion of entanglement, the 
main ingredient for quantum computation. It is entanglement makes the drastic dif- 
ference as compared with classical computation. The exponential dependence of the 
dimension on n is another reason for the attractiveness of quantum computing. The 
underlying picture in this case is therefore given by a multi-particle quantum theory, 
recall for example the discussion of the Einstein- Podolsky-Rosen paradox. A reconcilia- 
tion of our discussion with the tensor product structure will need further investigations 
but again, our main motivation was to propose a possible physical realization of unitary 
operators associated to elementary quantum gates. 

Our approach may also be adapted to the situation where we consider nonrelativistic 
electrons with spin moving through the wires and with boundary conditions possibly 
allowing for spin flips. Now for n wires the Hilbert space of incoming plane waves 
at fixed energy will have dimension equal to 2n. As mentioned in the introduction 
we will not discuss the experimental feasibility of designing such vertices with a given 
boundary condition. Also we will not discuss to what extent the dynamics given by 
the Laplace operator A(A, B) may serve as an approximation to the dynamics of a real 
particle (like an electron, say) in concrete wires, where for example one is also confronted 
with the problem of localization and the resulting exponential decay in space of some 
eigenfunctions due to the presence of impurities (see e.g. [^5[ for a discussion on this 
point) and which may prevent the transmission (see [pi]]). 

Example 1. In what follows the energy E will he fixed at Eq. Consider the one- 
dimensional unitary matrix given as a phase exp(ix) (n = 1). Then the Robin boundary 
condition cos(p ^(0) + sine/; ^'(0) = with (f) chosen such that 

,v cos <j) — iy/Eo sin 6 

e^ = ■^= 

cos (p + iV^o sin <j) 

solves this problem for one wire. Next consider the 2x2 unitary, idempotent so-called 
Hadamard matrix {n = 2) 

-(' ' 

V2\l -1 



Then the boundary conditions are given by 



A 



B 




Note that both A and B have vanishing determinant. 

Finally consider the real, unitary 4x4 matrix (the XOR gate, which is a special 
Toffoli gate, also called a controlled- NOT ("CNOT")) {n = 4) 



/I 









1 

0/ 



The corresponding boundary conditions are now given as 



A 



B 



/ \ 

i -i 



V 



/ 



1 



2i^/Ei 







/ 2 \ 

2 

11 

V 1 1 / 



Again both A and B have vanishing determinant. 

The S-matrix at a given energy and resulting from connecting such (elementary) 
gates is obtained from the S-matrices at the same energy of these individual gates and 
the lengths connecting them by what we cahed the generalized star product [po| |23|| , 
because it generahzes the star product in classical network theory (see p8| , |39|]). Al- 
though associative this composition rule is nonlinear due to the presence of reflection 
amplitudes. In the context of potential (quantum) scattering theory on the line an 
equivalent formulation of the star product also figures under the name Aktosun formula 
(see e.g. jT], ^, 43, 44, |l^) and there it is an easy consequence of the multiplicative 
property of the fundamental solution of the Schrodinger equation for a given energy and 
which is a 2 X 2 matrix (see [^] ) . In the theory of mesoscopic systems and multichannel 
conductors this has also been known for a long time and there it is called the multiplica- 
tivity of the transfer matrix, which is conjugate to the fundamental solution (see e.g. 
[53, 11, 12, 13, 14, 29, 49, ^, |l^)). In this context we would also like to mention that 
the problem of how to connect gates in a concrete way without feedback problems (back 
propagation) due to the presence of refiection amplitudes seems to be well known (see 
e.g. y, p. 3458], [^]) but to the best of our knowledge as yet no concrete proposal has 
been made, even on the conceptual level. Thus so far ordinary matrix multiplication has 
been used as the composition rule when connecting several gates (see e.g. Q and the 
references quoted there) . This is in analogy to the notion of straightline programming 
in information theory. 



3 Experimental determination of the S-matrix 

In this section we will address the general problem to what extent scattering experiments 
fix the S-matrix. In what follows the quantum wire will be arbitrary, i.e. it may contain 
an arbitrary number of external and internal lines and with an arbitrary number of 



vertices (see [20| for details). First we recall that the matrix elements Sjk{E) introduced 
in (|^) have the physical interpretation of transmission coefficients {j 7^ k) and reflection 
coefficients (j = k) into channel j for an incoming plane wave of energy E in channel 
k. Hence their absolute values are observable. On the other hand if one prepares an 
incoming wave packet (j){-; E) to be of the form 



ix;E) = AfcC" 



/Ex 



with arbitrary complex amplitudes AfcGC(l</c<n), then the resulting transmission 
amplitudes into channel j are of the form 



Y^XkSjkiE), (8) 



fc=i 



whose absolute values in principle are observable. But an easy argument implies the 
following. By a relabeling of the channels if necessary the S-matrix may be decomposed 
into a direct sum as 

Sa,b{E) = (B^^^S^siE) (9) 

for all -E > 0, where M is supposed to be maximal with 1 < M < n. Here the S^^{E) 
are unitary km x km matrices with "^km = n. Indeed, since Sa,b{E) is real analytic 
m. E > Q any of its matrix elements is either identically equal to zero or vanishes on at 
most a denumerable set without accumulation points. This implies that such a maximal 
decomposition (P) exists and is unique. From the point of view of boundary conditions 
this means that {A, B) may be chosen to have a corresponding maximal decomposition. 
This means that actually one has a disconnected graph (see |^0| for further details). 
Now experiments fix Sa,b{E) up to a phase factor matrix exp (ix{E)), where 

exp (ixXE) = e*'=i exp iiXm{E))I. (10) 

So assume now that experimentally one has measured the S-matrix for all energies in 
the form S{E) which agrees with Sa,b{E) up to a phase factor matrix exp{ix{E)), i.e. 

the relation Sa,b{E) = exp {ix{E)) S (E) holds for as yet undetermined phase factors 
exp {ixm{E)), 1 < m < M. In particular note that the decomposition (y) is observable, 
i.e. S{E) has a decomposition of the same form. Since we know that Sa,b{E) is real 
analytic in \/E > we may as well assume S{E) to be real analytic by fitting with a 
phase factor matrix if necessary. But then the phase factor matrix exp(ix(-E')) is also 
real analytic. In particular if the graph has one vertex only then by (0) this results in 
the following relation for the phase factor matrix 



exp {ixiE)) = S{E)-' (Ve - y%) exp {ix{Eo))S{Eo) + {Ve + y%) ^ 

■ (^{Ve + y%) exp {ix{Eo))S{Eo) + {Ve - y%)) (11) 

for all E and Eq. Now fix Eq. If there is a solution, then (^) shows that exp {ix{E)) is 
completely determined in terms of the "initial condition" exp (ix(-E'o)) for all E. On the 



other hand in order to find a solution one has to choose exp (ix(-Eo)) in such a way that 
the r.h.s. of (pl|) is a phase factor matrix for all i? > in the sense of (0). However, 
this procedure does not lead to a unique solution in general as may be seen from looking 
at Robin boundary conditions given as 

Ajk = ^jk cos 0fc, Bjk = djk sin (/)fc 

resulting in an S-matrix of the form 

cos(/>fc -i\/^sin(/>fc 

bjk{h) = -djk -— — -^=— — -. 

cos (pk + X V -C/ sm (pk 

Then one can essentially only achieve 

Sjk{E) = 5jk 

by experiments and any initial condition leads to a solution. This general lack of unique- 
ness is reminiscent of a similar situation in scattering theory ([^, ^, 26, 1^, 25], see also 
e.g. 



However, there is a way out to find the phases of the S-matrix. The idea is to enlarge 
the graph by components for which the associated S-matrices S'^"^(£') are supposed to 
be known. The enlarged graph will give rise to an S-matrix S^'^™{E), which is obtained 
from the original S-matrix S{E), S'^^^{E) and the lengths of the new internal lines 
(which also are supposed to be known) via the generalized star product [gO|, |2^. By 
varying S^^^{E) and measuring the absolute values of the matrix elements of S^^^{E), 
we may infer the phases of the matrix elements of S{E). This strategy is inspired by 
a familiar procedure in quantum computation, where one analyzes the effect of a black 
box (an "oracle" ) by combining the black box with another box, whose effect is known 



(see e.g. [34| and the references quoted there). 

We start with the case where we want to determine the phase of the reflection 
amplitude Su^E). If Sii{E) = 0, there is nothing to prove, so we will assume Su^E) / 0. 
On the external line labeled by i we introduce an additional vertex at a distance Oj from 
the vertex where the line ends. There we introduce arbitrary boundary conditions which 
result in an arbitrary 2x2 unitary S-matrix S^^'^{E) by our previous discussion. By the 
techniques of the generalized star product the resulting S-matrix S^°^{E) is obtained 
from S{E), S^^^{E) and the length Oj such that in particular 

SfriE) = U22 + U2iSii{E){l - Sii{E)Uii)-^Ui2- (12) 

Here the Uku ^) ^ = 1)2 are the matrix elements of the unitary 2x2 matrix 

In the article pC| ] using a "Born series" expansion we explain the physical intuition 
behind the generalized star product, from which in particular (11^ ) may easily be ob- 
tained. In |23] we shall provide a rigorous proof that the generalized star product indeed 



provides the right factorization. Since S^^^{E) and Oj are supposed to be known, U is 
also known and may be chosen arbitrary in the unitary group U(2). Thus Sff'^{E) is 
a function of Sui^E) and U and by measurements \Sf™{E)\ and \Sii{E)\ are known. 
Actually Sf™{E) = exp{2i^/Eai) Sii{E) for the special case when 



S^"^(^) = So 



1 

1 



the "unit matrix" with respect to the generaUzed star product. 
Write Sii{E) = exp(i(/.) \Sii{E)\. Thus 

\Sfr{E)\ = \U22 + U2i\Sii{E)\e'^{l - \Sii{E)\e'^Uii)-^Ui2\ 

is known for all U and the aim is to determine (/>. Now any [/ G U(2) may be written as 

/ e*(x+r)p -e^(x-)(i _ p2)i/2 X 

^ - \^ ei(x+«)(l_p2)i/2 e*(^-")p ) 

with < p < 1. Therefore 

p- (1 -p2)|5,,(ii;)|(e-(^+x+-) - \S,,{E)\pr^\ (13) 

is known. Now choose p so small that 2p^\Sii{E)\ < \Sii{E)\ — p. Then (|l^) is easily 
seen to be maximal for exp(i((;/) + x + r)) = — 1. Thus tuning x and r and hence S^'^^^E) 
while keeping p fixed, such that \Sf™{E)\ becomes maximal, fixes (p as was claimed. An 
alternative way to determine (p is to consider 

5f/-(i?) = S^j{E){l - Sii{E)Uii)-^Ui2 

with i ^ j under the same situation. This procedure therefore is only possible if the 
graph considered has at least two external lines. Now 



\Sf^'"{E)\ = \S,j{E)\{l - p'yi' e-H*+>^+-^ - \Su{E)\p 

becomes maximal when exp(i((/) + x + T)) = 1 and minimal when exp(i(</> + x + T)) = —1. 
We turn to a discussion of the phases of the transmission amplitudes Sij{E), i ^ j. 
As discussed above (see (^)), by a suitable preparation of the incoming state as a super- 
position of plane waves of energy E in the incoming channels i and j, we can measure 
\Sij{E) + A Sii{E)\ for all complex A. Since Sii{E) is has already been determined, this 
fixes Sij{E) for i / j provided Su^E) ^ 0. Now assume Su^E) = 0. Then we proceed 
as follows. By the same procedure as above we insert a vertex on the line labeled by 



i at a distance a[ with an S-matrix S^^^{E)' . In analogy to (12) we then obtain an 
S-matrix S{E)' with SniE)' = U22 = S22'^{Ey , which is non vanishing provided p' > 0. 
Then we can measure Sij{E)' with the arguments used above. Now S{E)'^,- converges 
to eyiY>{i^fEa[) Sij{E) when S'^^^{E)' converges to 5*0, so this determines Sij{E). The 
appearance of the phase factor exp(iv^ai) is due to the fact that the insertion of the 
new vertex shifts the origin of the line labeled by i by the distance a[. This concludes 
our proof that all matrix elements Sij{E) may be obtained from scattering experiments 
performed in this way. Note that for this procedure only scattering experiments at fixed 
energy are necessary. 

4 Description in terms of von Neumann's extension theory 

With the observations of the Section 2 we also may describe A{A, B) defined on a graph 
with a single vertex and the resulting S-matrix from the viewpoint of von Neumann's 
extension theory (see e.g. [^]). This extends the discussion in Appendix A of ||20|| . 
According to von Neumann's theorem any self-adjoint extension A of A*^ can be uniquely 
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parameterized by a linear isometric isomorphism W : Ker(— A*^ — i) — > Ker(— A'' + i) 
according to the formula 

V{A) = |'0 + V'+ + >V^+| i^eV{A°),i;+ GKer(-A°'^-i)|, 
-A(V' + tp+ + >VV'+) = -A°^ + #+ - i>V^+. 

We will describe W in terms of a matrix W by choosing particular bases in Ker(— A^ —i) 
and Ker(— A^ + i) respectively. 

Namely we choose Uj € Ker(— A'^ — i) and Vj G Ker(— A'^ + i), j = 1, . . . , n to be 
given as 



of 



(%(x)), = 5,fc2V4e73(-i+*)-, (^;^.(^)), = 5^.,2i/4e75(-i-*)^ fc = 1, . . . ,n. 

One can easily verify that {u}^^^ and {v}"^]^ are orthonormal bases for Ker(— A^' — i) 

and Ker(— A*^ + i), respectively. Thus W is the unitary matrix representation of W 
with respect to these bases {u}'j^i and {v}'j^i, i.e. 

n 

y^uj = ^WkjVk. 

k=l 



In Appendix A of |20| it was shown that for A = A(^, B) the matrix W = Wa,b was 
given as 

Wa,b = - (a - -^(1 + i)B^ (^A + -^(-1 + i)B^ . 

But this means that given any unitary matrix W describing a selfadjoint extension in 
the sense of von Neumann in the basis described above, and using the pair {A' , B') given 
as {A', B') = {CA, CB) with C = V2{A - 1/^/2(1 + i)B)-^ = ^/2{A - exp(i7r/4)5)-i to 
describe the boundary conditions, we have W = —1 — iB' as well as W = 1— \f2N + B' 
giving A! and B' in terms of W as 

A! = -e-'^'I^W ^e^'lH 
B' = i{W + I). 

Therefore the S-matrix Sy^{E) for a selfadjoint extension given in terms of W, and 
which has a matrix representation W in the above bases, takes the form 



■ ((-e-^"/^ + ^/E)W + (e*"/^ + /E)i) . (14) 

As a consistency check note that Sy^;{E) is unitary for all £J > since W is unitary. 
Also W = —I corresponds to Dirichlet boundary conditions with Sa=i,b=o{E) = —I for 
all i? > 0. This last relation may of course be inverted to give the matrix W in terms 
of the S-matrix at any energy. Indeed 



W = ((e*"/^ - \/;E)5(^) + (e^"/^ + /E)i) 

• ((e"^"/^ + Ve)S{E) + {e-'^l^ - ^fE)\\^ '^ . (15) 
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In particular the right hand side is independent of the energy E so this may serve 
as another test that the S-matrix indeed results from boundary conditions in the way 
described above. 

Acknowledgements: The authors would like to thank David Tomanek for infor- 
mative discussions concerning nanotubes and an anonymous referee for helpful critical 
remarks. 
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